THE UNIVERSE EVOLUTION IN EXPONENTIAL GRAVITY 
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A generic feature of viable exponential _F(i?)-gravity is investigated. An additional modification to stabilize the effective dark 
energy oscillations during matter era is proposed and applied to two viable models. An analysis on the future evolution of the 
universe is performed. Furthermore, a unified model for early and late-time acceleration is proposed and studied. 



1 Introduction 

The current acceleration of our universe is supported 
by various observations and cosmological data pQ. 
There exist several descriptions of this acceleration 
and, among them, the most important are given by 
the ACDM Model of General Relativity (GR) plus Cos- 
mological Constant and by the modified gravitational 
theories, like F(i?)-gravity [2J. Here, we study some as- 
pects and generic feature of viable F(i?)-gravity models 
which reproduce the physics of ACDM Model, namely 
power-law and exponential gravity. We show that the 
behavior of higher derivatives of the Hubble parame- 
ter may be influenced by large frequency oscillations of 
effective dark energy during matter era, which makes 
solutions singular and unphysical at a high redshift. 
As a consequence, we examine an additional correction 
term to the models in order to remove any instability 
with keeping the viability properties. Some comments 
about the growth of matter perturbations of corrected 
models are given and future universe evolution is ana- 
lyzed. We also propose a way to describe the early time 
acceleration of universe by applying exponential grav- 
ity to inflationary cosmology and a unified description 
between inflation and current acceleration is presented. 
This work is based on Refs. [21 H]. We denote the grav- 
itational constant 87tGn by n 2 . 

2 Realistic F(i?)-models in the FLRW universe 

The action describing F(i?)-gravity is given by 

F(R) i ^(matter 



2n 2 



(1) 



where F(R) is a generic function of the Ricci scalar R 
only, g is the determinant of the metric tensor g^ v and 
^(matter) j g fa e mat ter La grangian. In realistic models 



reproducing the standard cosmology of GR with a suit- 
able correction to realize current acceleration and/or 
inflation one has F(R) = R + f(R) and the modifica- 
tion of gravity is encoded in the function f(R). In this 
case, we may use an effective fluid-like representation 
of modified gravity and in FLRW space-time described 
by the metric ds 2 = —dt 2 + a(t) 2 dx 2 , the equations of 
Motion (EOM) assume the form 



Pde + Pm 
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(2) 

Here, p m and P m are the energy density and pres- 
sure of standard matter and pde and Pde are the 
effective dark energy density and pressure given by 
modified gravity. In order to study the dynamics of 
F(R) gravity models, we may introduce the variable 
Vh(z) = /9DE/Vm(o) as a function of the red shift z, 
where p m (o) = 3m 2 /« 2 is the matter energy density 
at the present time, being rn 2 a mass scale. From the 
EOM one derives an equation on the type [5] 

y'^(z) + Ji(z)y' H {z) + J 2 (z)y H {z) + J 3 (z) = 0, (3) 

where J\{z), J 2 (z) and Jz(z) depend on the F(R)- 
model. In what follows, we will consider viable models 
representing a realistic scenario to account for dark en- 
ergy, in particular a suitable reparameterization of the 
Hu-Sawicki model [6], 



F(R) = R-2Ail- 
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OR/A) 4 + 1 

and a simple example of exponential gravity [7] , 



F(R) =R-2A 1 - e 



-R/A 



(4) 



(5) 



This models are very carefully constructed such that in 
the high curvature regime F(R ^> A) = R — 2A. Here, 
A plays the role of the Cosmological Constant which 
decreases in the flat limit (F(R — > 0) — > 0) where the 
Minkowski solution is recovered. 
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Figure 1 : Evolution of the jerk parameter in pure exponential model (a) and in exponential model with correction 
term (b) overlapped with the ones in ACDM Model. 



3 DE-oscillations in the matter dominated era 

Since in matter dominated era R — 2>fh 2 (z + l) 3 , one 
may locally solve Eq. (|3|) around z+Sz, where \Sz\ <C z. 
The solution reads 

y H (z + Sz) = a + b Sz + C e W (-odV^)* _ 

(6) 

where Cq is a generic constant and arj, bo, a and 
/3 are constants depending on the model. It turns 
out that the dark energy perturbations remain small 
around z as soon as the exponential term does not 
diverge in expanding universe, when Sz < 0. The ex- 
plicit conditions (1-F'(R))F"'(R)/(2F"(R) 2 ) > -7/2 
and 1/(RF"(R)) > 12 are found for stable matter 
era [3]. In the case of realistic gravity described by 
the models (HI)-©, where a = A/(3m 2 ), it is easy to 
see that such conditions are well satisfied. Further- 
more, since in both of the models F"(R) — > + , we 
understand that the discriminant of Eq. © is nega- 
tive and dark energy oscillates with frequency J-(z) ~ 
(R * F (R)) 1 /(z + 1). Since by increasing the cur- 
vature F"{R) tends to zero, for large values of the red- 
shift the dark energy density oscillates with a high fre- 
quency and also its derivatives become large and may 
be dominant in some derivatives of the Hubble parame- 
ter which approach an effective singularity making un- 
physical the solution. Moreover, it may be stated that 
the closer the model is to the ACDM Model (i.e., as 
much F"(R) is close to zero), the bigger the oscillation 
frequency of dark energy becomes. As a consequence, 
despite the fact that the dynamics of the universe de- 
pends on the matter and the dark energy density re- 
mains very small, some divergences in the derivatives 
of the Hubble parameter can occur, showing a different 
feature of this kind of models with respect to the case of 
GR with an undynamic Cosmological Constant. Since 

lr The form of the jerk parameter as a function of time is j(t) 
2 We have used Mathematica 7 ©. 



in the models (HJ)-© the approaches to the Cosmologi- 
cal Constant are different from each other, we may say 
that dark energy oscillations in matter era are a generic 
feature of realistic F(R) gravity, in which the cosmo- 
logical evolution is similar to those in ACDM Model. 
In order to remove the divergences in the derivatives 
of the Hubble parameter, we can introduce a function 
g(R) for which the oscillation frequency of the dark en- 
ergy density acquires a constant value J~(z) = l/VS, 
S > 0, during matter era, stabilizing the oscillations 
with the use of a correction term. This term has been 
derived as 

/ R \i/3 
^) = -7a(^) , (7) 

where 7 = 1.702 5, so that (R * g"(R)y 1/2 /(z + 1) 
when R = 3fh 2 (z + l) 3 . We explore the models 
(HI)-© with adding these correction. The effects of 
the compensating term vanish in the de Sitter epoch, 
when R — AL and the models resemble to a model 
with an effective cosmological constant, provided that 
7 <C (to 2 /A) 1 / 3 . For example, a reasonable choice is 
to put 7 = 1/1000, which corresponds to the oscilla- 
tions frequency of the dark energy during matter era 
T = 41.255 with a period T ~ 0.152. In order to 
appreciate the role of the proposed correction, let us 
consider the jerk parameter j(z), which depends on 
the second derivative of i?(^lj, for the exponential 
model ©, whose numerical extrapolation^ is plotted 
and overlapped with its behavior in ACDM Model on 
the left side of Fig. Q] Despite the fact that the dark 
energy density pde — A/k 2 is negligible at high red 
shift, the effects of its oscillations become relevant in 
the jerk parameter, which grows up with an oscillatory 
behavior and diverges at z ~ 2.8. On the right side of 
the figure, we plot the graphic of the jerk parameter 
for the same model with the addition of the correction 

a'" (t) I (a(t)H(t)), a(t) being the scale factor of the universe. 
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term g{R) previously discussed. Also in this case the 
jerk parameter oscillates with the same frequency of 
the dark energy, but here such oscillations have a con- 
stant frequency and do not diverge at high red shift. 
It is also important to stress that the term added 
to stabilize the dark energy oscillations in the matter 
dominated epoch does not cause any problem to the 
good proprieties of the models ((I])-©, which continue 
to satisfy all the cosmological and local gravity con- 
straints. It has been demonstrated in several papers 
that in order to evade the solar system tests [8] or 
avoid large corrections to Newton Law at large cosmo- 
logical scales [9] , a modified gravity theory must satisfy 
the following relation 

3k(H-«)»>. 

when R assumes the typical curvature value of the 
constant background around a matter source. For 
the models ((4j)- (|5]) without correction terms, this con- 
dition is always well satisfied, due to the fact that 
F"(R > A) ~ + . Next, if we consider the addition of 
g(R), whose second derivative is dominant, the left side 

of the letter equation reads 3 4 / 3 m 2 / (2Aj) (R/m 2 ) 2/3 . 
Despite the fact that this quantity is smaller than the 
one obtained for the pure models ©-©, it still re- 
mains sufficiently large and the correction to the New- 
ton Law is very small. For example, at the typical 
value of the curvature in the solar system, namely 
R ~ 10~ 61 eV 2 , it corresponds to 2 x 10 6 , which is 
very large effectively Concerning the matter instabil- 
ity [10], this might also occur when the curvature is 
rather large, as on a planet (R ~ 10 _38 eV 2 ), as com- 
pared with the average curvature of the universe today 
(R ~ 10 _66 eV 2 ). This instability is avoided as soon 
as 1/(RF"(R)) is much larger than one on the planet 
surfaces. For our correction term g(R), this quantity 
evaluated at R ~ 10 _38 eV 2 corresponds to 4 x 10 21 and 
also in this case we do not have any particular problem. 

4 Growth index 

After the discovery of cosmic acceleration, a lot of the- 
ories for the dark energy have been proposed, and it 
has become very important to find a way to discrimi- 
nating among these theories, most of them exhibit the 
same expansion history. In this context, the charac- 
terization of growth of the matter density perturba- 
tions is very significant in modified gravity, since mod- 
els with the same cosmological evolution have a dif- 
ferent cosmic growth history. In order to execute it, 
the so-called growth index j(z) is useful [UJ. Un- 
der the subhorizon approximation, the matter den- 



sity perturbation S = 5p m /p m satisfies the equation 
6 + 2H8 — 47rG c ff (a, k)p m S — with k being the co- 
moving wavenumber and G e ff(a, k) being the effective 
gravitational "constant" depending on the model and 
expressed as a function of k and the scale factor a, 
such that k 2 /a 2 ^> H. The equation for matter per- 
turbation has an equivalent description in terms of the 
growth rate / g = dlnS/dhia. The growth index j(z) 
is defined as is defined as the quantity satisfying the 
relation f g {z) = fi m (» 7(z) with Q m (z) = £,TrGp m /3H 2 
being the matter density parameter. The growth index 
cannot be observed directly, but it can be determined 
from the observational data of both the growth factor 
and the matter density parameter at the same redshift 
and this is the reason for which in the last period its 
calculation in different theories has become of a great 
interest. In Ref. [4] various parameterizations for the 
growth index 7(2) have been studied for the models 
©-© with correction term g(R). It has been found 
that the Ansatz which better fits the (numerical) so- 
lution of the growth rate for the considered models is 
given by 7 = 70 + 712, where 70,1 are constants. 

5 Future universe evolution 

The confrontation of the models ©-© with SNIa, 
BAO, CMB radiation and gravitational lensing has 
been executed in the past in several works and all 
the numerical analysis have demonstrated that they 
are consistent with the last observational data com- 
ing from our universe. In particular, they repro- 
duce the correct amount of dark energy today J^de = 
0.721±0.015 [SI [JJ]. Moreover, the addition of cor- 
rection term © does not have any influence at the 
present day range of detection. Starting from Eq. ©, 
we are able to study perturbations around the de Sit- 
ter solution R d § = 4A of our universe today in the 
considered models. Performing the variation with re- 
spect to Uh(z) = Ho + Vi(z) with yo = R d $ / 12m 2 and 
I I <C 1 and assuming the contributions of radia- 
tion and matter to be much smaller than yo, one finds 
the following solution for yi(z), 

Vl (z) = C (z + l)K3±v^) + A A {z + 1)3 , (9) 

where Co is a generic constant and C = 1 + (1 — 
F'(R))/{RF"(R)) and /3 = —4 + 4F' (R) / (RF" (R)) , 
the both evaluated at R = R d s- The de 
Sitter solution is stable provided by condition 
F'(Rds)/{{Rds)F"{R dS )) > 1, which is a well know 
result. It has also been demonstrated that since in real- 
istic F(i?)-gravity models for the de Sitter universe < 
Rd S F"(R ds ) « 1, one has F'(R ds )/(R ds F"(R ds )) > 
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25/16, giving the negative discriminant of Eq. © and 
an oscillatory behavior to the dark energy density dur- 
ing this phase. Thus, in this case the dark energy 
Equation of State (EoS) parameter lode = Pde/pde 
oscillates infinitely often around the line of the phan- 
tom divide wde = — 1 j!3) . These models possess one 
crossing in the recent past, after the end of the matter 
dominated era, and infinite crossings in the future, but 
the amplitude of such crossings decreases as (z + l) 3 / 2 
and it does not cause any serious problem to the accu- 
racy of the cosmological evolution during the de Sitter 
epoch which is in general the final attractor of the sys- 
tem. However, since the existence of a phantom phase 
can give some undesirable effects such as the possibility 
to have the Big Rip [T4] or Big Rip scenario with the 
disintegration of bound structures [15] , it is important 
to observe that in the models (|4j)- (J5j) with the adding 
modification (0, the effective EoS parameter of the 
universe 
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PDE + Pn 
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2(z + l)dH{z) 
3H{z) dz~ 



(10) 



never crosses the phantom divide line in the past, and 
that only in the very far future, when z is very close to 
— 1, it coincides with lode and the crossings occur. Let 
us consider for example the exponential case ([5]) plus 
term Q. The numerical evaluation of Eq. © leads to 
Hubble parameter 



H(z) = y/m* [y H (z) + (z + m 



(11) 



and therefore u> e g(z). We depict the corresponding cos- 
mological evolution of Lo e g(z) for — 1 < z < 2 in Fig. 2. 
We can see that u> e g(z) starts from zero in the matter 
dominated era and asymptotically approaches -1 with- 
out any appreciable deviation. Only when z ~ —0.90, 
u> e g(z) crosses the line of phantom divide, how it is 
shown on the right. By considering the scale factor 
a(t) = exp(H t), where H ~ 6.3 x 10- 34 eV _1 is the 
Hubble parameter of the de Sitter universe, we may 
conclude that the crossings appear at 10 26 years. 




6 Exponential gravity for early and late-time 
cosmic acceleration 

One of the most important goals for modified grav- 
ity theories is that all the processes in the expan- 
sion history of the universe are realized, from infla- 
tion to the late cosmic acceleration epoch, without in- 
voking the presence of dark components. Models of 
the type ((SJ) may be combined in a natural way to 
obtain the phenomenological description of the infla- 
tionary epoch, mimicking a smooth version of F(R) = 
R-2A [l - e- R / A ] -Aj O(R-Ri), where R t is the tran- 
sition scalar curvature for which a suitable cosmological 
constant Aj switches on producing inflation. Following 
the first proposal of Ref. [TB] , we consider the form of 
F(R) with a natural possibility of a unified description 
of our universe 



F(R) 



R-2A11 



+7 



A, 



1 



e-(*)' 



R? 



R a 



(12) 



where i?j and A; are the typical values of transition 
curvature and expected cosmological constant during 
inflation, respectively, and n is a natural number larger 
than unity (here, we do not write the correction term 
for the stability of oscillations in the matter dominated 
era). In Eq. (|T2j, the last term j(l/R"' 1 )R a 1 where 
7 is a positive dimensional constant and a is a real 
number, works at the inflation scale i£; and is actu- 
ally necessary in order to realize an exit from infla- 
tion. An accurate analysis of this model shows that if 
a > 1 and n > 1, we avoid the effects of inflation at 
small curvatures and it does not influence the stabil- 
ity of the matter dominated era. The de Sitter point 
i?ds which describes the early time acceleration exists 
if Ri = i?as and it reads as i?ds = 2/(7(2 — a) + 1) 
under the condition (R d s/Ri) n 3> 1. Thus, the infla- 
tion is unstable if aj(a — 2) > 1. We may evaluate 
the characteristic number of e-folds during inflation 
N = log[(zj + l)/(z c + 1)], where z\ and z c are the 
rcdshifts at the beginning and at the end of early time 
cosmic acceleration. Given a small cosmological per- 
turbation y\{zi) at the redshift z- x to the effective dark 
energy density of inflation y — i?ds/(12m 2 ) , so that 
y(z) = ya + yi{z), we have from Eq. ^ avoiding the 
matter contribution 



yi(zi) = C {z { + i) x , 



Figure 2: Effective EoS parameter for — 1 < z < 2 in wrtn 
exponential model. On the right, the region — 1 < z < 
—0.9 is stressed. 




l6F'(R dS ) 
Rd S F"(R ds ) 
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(14) 
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Here, x < if the de Sitter point is unstable. As a 
consequence, the perturbation yi(z) grows up in ex- 
panding universe as yi(z) = yi(zi) \{z + l)/(zi + ^)) x ■ 
When j/i (z) is on the same order of the effective mod- 
ified gravity energy density yo(z) of the de Sitter solu- 
tion describing inflation, the model exits from inflation. 
In this way, we can estimate the number of e-folds dur- 
ing inflation as 



N 



1 



loe 



2/0 



(15) 



A value demanded in most inflationary scenarios is 
at least N = 50-60. A classical perturbation on the 
(vacuum) de Sitter solution may be given by the pres- 
ence of ultrarelativistic matter in the early universe. 
Thus, we can predict the numbers of e-folds for the 
model (fT2|) by computing x in Eq. (Ti"4l) , which depends 
on model parameters. Let us consider an example 
by supposing that ultrarelativistic matter/radiation is 
10 6 times smaller than that of effective modified grav- 
ity energy density (~ Ai) at the beginning of infla- 
tion. By choosing the viable model parameters n — 4, 
Ri = 2Aj, 7 = 1 and a — 8/3, such that the curva- 
ture at the inflation results to be i?ds = 6Ai, one has 
x = —0.218 and, by using Eq. ([15]). we obtain N ~ 64 
and z c = e~ 64 (z; + 1) — 1. The numerical extrapo- 
lation of Eq. (|3) gives us ynizc) — 0.88j/ff(zi) and 
R(z e ) = 0.853i?dS: confirming that the effective modi- 
fied gravity energy density and the curvature decrease 
at the end of inflation. After that, it is reasonable 
to expect that the physical processes described by the 
ACDM Model are reproduced by the first term of mod- 
ification to gravity. 

7 Conclusions 

In this work, we have analyzed some feature of realistic 
i 7 '(i?)-gravity mimicing the physics of ACDM Model. 
Two viable models have been considered. We have 
seen that since the corrections to GR at small curva- 
tures may lead to undesired effects at high curvatures, 
some adding modifications are required. In particular, 
we have reconstructed a compensating term in order 
to stabilize the oscillations of the effective dark energy 
at high red shifts with retaining the viability propri- 
eties. Furthermore, we have briefly discussed about the 
growth index and an investigation on the cosmological 
evolution of the universe has been executed showing 
that the effective crossing of the phantom divide, which 
characterizes the de Sitter epoch, occurs in the very 
far future. We also have considered the inflationary 
cosmology in exponential gravity. It has been derived 
the correlation between the e-folds during inflation and 
the model parameters. Since at the end of inflation the 



effective modified gravity energy density and the cur- 
vature decrease, we may conclude that it is possible to 
acquire a gravitational alternative scenario for a uni- 
fied description of inflation in the early universe with 
the late-time cosmic acceleration. 
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